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ON . ABSTRACT 

ON. 

-*— » , The energy density of a scalar field (j) with potential V(4>) oc 4>~ a , a > 0, behaves like 

/•""n j a time-variable cosmological constant that could contribute significantly to the present 

energy density. Predictions of this spatially-flat model are compared to recent Type la 

fSj ' supernovae apparent magnitude versus redshift data. A large region of model parameter 

space is consistent with current observations. (These constraints are based on the exact 
J> ■ scalar field model equations of motion, not on the widely used time-independent equa- 

tion of state fluid approximation equations of motion.) We examine the consequences of 

l/~j , also incorporating constraints from recent measurements of the Hubble parameter and 

the age of the universe in the constant and time- variable cosmological constant models. 

q^ ' We also study the effect of using a non- informative prior for the density parameter. 

Ov 

Subject headings: supernovae: general — cosmology: observations — large-scale structure 
*j"H of the universe 

O' 

a: 

f~ ' 1. Introduction 

>: 

Current observations are more consistent with low-density cold dark matter (CDM) dominated 
cosmogonies. For recent discussions see, e.g., Ratra et al. (1999a), Kravtsov & Klypin (1999), 
Doroshkevich et al. (1999), Colley et al. (2000), Freudling et al. (1999), Sahni & Starobinsky 
(1999), Bahcall et al. (1999), and Donahue & Voit (1999). The simplest low-density CDM models 
have either flat spatial hypersurfaces and a constant or time-variable cosmological "constant" A 
(see, e.g., Peebles 1984; Peebles & Ratra 1988, hereafter PR; Efstathiou, Sutherland, & Maddox 
1990; Stompor, Gorski, & Banday 1995; Caldwell, Dave, & Steinhardt 1998), or open spatial 
hypersurfaces and no A (see, e.g., Gott 1982, 1997; Ratra & Peebles 1994, 1995; Gorski et al. 
1998). 

While these models are consistent with most recent observations, there are notable exceptions. 
For instance, recent applications of the apparent magnitude versus redshift test based on Type la 
supernovae favor a non-zero A (see, e.g., Riess et al. 1998, hereafter R98; Perlmutter et al. 1999a, 
hereafter P99), although not with great statistical significance (Drell, Loredo, &i Wasserman 1999). 

On the other hand, the open model is favored by: 



• 



Measurements of the Hubble parameter Hq (= lOO/i km s 1 Mpc ) which indicate h = 
0.65 ± 0.1 at 2 a (see, e.g., Suntzeff et al. 1999; Biggs et al. 1999; Madore et al. 1999), and 
measurements of the age of the universe which indicate to = 12 ± 2.5 Gyr at 2 a (see, e.g., 
Reid 1997; Gratton et al. 1997; Chaboyer et al. 1998). This is because the resulting central 
-ffo^o value is consistent with a low-density open model with nonrelativistic-matter density 
parameter Qq ~ 0.35 and requires a rather large Qq ~ 0.6 in the flat-constant-A case. 

Analyses of the rate of gravitational lensing of quasars and radio sources by foreground 
galaxies which require a large f^o > 0.38 at 2 a in the flat-constant-A model (see, e.g., Falco, 
Kochanek, & Muhoz 1998). 

Analyses of the number of large arcs formed by strong gravitational lensing by clusters (Bartel- 
mann et al. 1998; Meneghetti et al. 1999). 

• The need for mild anti-biasing to accommodate the excessive intermediate- and small-scale 
power predicted in the DMR normalized flat-constant-A CDM model with a scale-invariant 
spectrum (see, e.g., Cole et al. 1997). 

While the time-variable A model has not yet been studied to the same extent as the open and 
flat- A cases, it is likely that it can be reconciled with most of these observations (see, e.g., PR; 
Ratra & Quillen 1992, hereafter RQ; Frieman & Waga 1998; Perlmutter, Turner, & White 1999b; 
Wang et al. 1999; Efstathiou 1999). 

We emphasize that most of these observational indications are tentative and certainly not 
definitive. This is particularly true for constraints derived from a x 2 comparison between model 
predictions and cosmic microwave background (CMB) anisotropy measurements (see, e.g., Ganga, 
Ratra, & Sugiyama 1996; Lineweaver & Barbosa 1998; Baker et al. 1999; Rocha 1999), see dis- 
cussions in Bond, Jaffe, & Knox (1998) and Ratra et al. (1999b). More reliable constraints follow 
from models-based maximum likelihood analyses of the CMB anisotropy data (see, e.g., Gorski et 
al. 1995; Yamamoto & Bunn 1996; Ganga et al. 1998; Ratra et al. 1999b; Rocha et al. 1999), 
which make use of all the information in the CMB anisotropy data and are based on fewer approx- 
imations. However, this technique has not yet been used to analyze enough data sets to provide 
robust statistical constraints. Kamionkowski & Kosowsky (1999) review what might be expected 
from the CMB anisotropy in the near future. 

In this paper we examine constraints on a constant and time- variable A that follow from recent 
Type la supernovae apparent magnitude versus redshift data and recent measurements of Hq and to- 
We focus here on the favored scalar field {4>) model for a time- variable A (PR; Ratra & Peebles 1988, 
hereafter RP), in which the scalar field potential V{<j)) oc (j)~ a , a > 0, at low redshift 1 . This scalar 



Other potentials have been considered, e.g., an exponential potential (see, e.g., Lucchin & Matarrese 1985; RP; 
Ratra 1992; Wetterich 1995; Fcrreira & Joyce 1998), or one that gives rise to an ultra-light pseudo-Nambu-Goldstone 
boson (see, e.g., Frieman et al. 1995; Frieman & Waga 1998), but such models are either inconsistent with observations 



field could have played the role of the inflaton at much higher redshift, with the potential V((f>) 
dropping to a non-zero value at the end of inflation and then decaying more slowly with increasing 
4> (PR; RP). See Peebles k Vilenkin (1999), Perrotta k Baccigalupi (1999) and Giovannini (1999) 
for a specific model and observational consequences of this scenario. A potential oc 4>~ a could arise 
in a number of high energy particle physics models, see, e.g., Binetruy (1999), Kim (1999), Barr 
(1999), Choi (1999), Banks, Dine, k Nelson (1999), Brax k Martin (1999), Masiero, Pietroni, k 
Rosati (1999), and Bento k Bertolami (1999) for specific examples. It is conceivable that such a 
setting might provide an explanation for the needed form of the potential, as well as for the needed 
very weak coupling of (f> to other fields (RP; Carroll 1998; Kolda k Lyth 1999, but see Periwal 1999 
and Garriga, Livio, k Vilenkin 1999 for other possible explanations for the needed present value of 
A). 

A scalar field is mathematically equivalent to a fluid with a time-dependent speed of sound 
(Ratra 1991). This equivalence may be used to show that a scalar field with potential V(4>) oc 4>~ a , 
a > 0, acts like a fluid with negative pressure and that the <fi energy density behaves like a 
cosmological constant that decreases with time. This energy density could come to dominate at 
low redshift and thus help reconcile low dynamical estimates of the mean mass density with a 
spatially-flat cosmological model. Alternate mechanisms that also rely on negative pressure to 
achieve this result have been discussed (see, e.g., Ozer k Taha 1986; Freese et al. 1987; Turner 
k White 1997; Chiba, Sugiyama, k Nakamura 1997; Ozer 1999; Waga k Miceli 1999; Overduin 
1999; Bucher k Spergel 1999). However, these mechanisms are either inconsistent or do not share a 
very appealing feature of the scalar field models. For some of the scalar field potentials mentioned 
above, the solution of the equations of motion is an attractive time-dependent fixed point (RP; 
PR; Wetterich 1995; Ferreira k Joyce 1998; Copeland, Liddle, k Wands 1998; Zlatev, Wang, k 
Steinhardt 1999; Liddle k Scherrer 1999; Santiago k Silbergleit 1998; Uzan 1999). This means 
that for a wide range of initial conditions the scalar field <f> evolves in a manner that ensures that 
the cosmological constant contributes a reasonable energy density at low redshift. Of course, this 
does not resolve the (quantum-mechanical) cosmological constant problem. 

There have been many studies of the constraints placed on a time-variable A from Type la 
supernovae apparent magnitude versus redshift data (see, e.g., Turner k White 1997; Frieman 
k Waga 1998; Garnavich et al. 1998; Hu et al. 1999; Cooray 1999; Waga k Miceli 1999; P99; 
Perlmutter et al. 1999b; Wang et al. 1999; Efstathiou 1999). However, except for the analysis of 
Frieman k Waga (1998), who use the earlier Perlmutter et al. (1997) data, these analyses have 
mostly made use of the time-independent equation of state fluid approximation to the scalar field 
model for a time-variable A (Perlmutter et al. 1999b and Efstathiou 1999 do go beyond the time- 
independent equation of state approximation by also approximating the time dependence of the 
equation of state, however, they do not analyze the exact scalar field model) 2 . A major purpose 



or do not share the more promising features of the inverse power law potential model. 

2 A similar criticism holds for most analyses of the constraints on a time-variable A from gravitational lensing 



of this paper is to use the newer supernovae data of R98 and P99 to derive constraints on a time- 
variable A without making use of the time-independent equation of state fluid approximation to 
the scalar field model. 

In the analyses here we use the most recent R98 and P99 data to place constraints on a constant 
and time-variable A. We note, however, that this is a young field and insufficient understanding 
of a number of astrophysical processes and effects (the mechanism responsible for the supernova, 
evolution, environmental effects, intergalactic dust, etc.) could render this a premature undertaking 
(see, e.g., Hoflich, Wheeler, & Thielemann 1998; Aguirre 1999; Drell et al. 1999; Umeda et al. 
1999; Riess et al. 1999; Wang 1999). Furthermore, other cosmological explanations (time- variable 
gravitational or fine structure "constants") could be consistent with the data (see, e.g., Amendola, 
Corasaniti, & Occhionero 1999; Barrow & Magueijo 1999; Garcia-Berro et al. 1999). 

In addition to analyzing just the supernovae data, we also incorporate constraints from recent 
measurements of Hq and to an d derive a combined likelihood function which we use to constrain both 
a constant and a time-variable A. We also examine the effect on the model-parameter constraints 
caused by varying the prior for £Iq. 

We emphasize that the tests examined in this paper are not sensitive to the spectrum of inho- 
mogeneities in the models considered. This spectrum (possibly generated by quantum mechanical 
zero-point fluctuations during inflation, see, e.g., Fischler, Ratra, &: Susskind 1985) is relevant for 
some of the other tests mentioned above (CMB anisotropy, anti-biasing, etc.). 

In §2 we summarize the computations. Results are presented and discussed in §3. Conclusions 
are given in §4. 



2. Computation 

We examine three sets of supernovae apparent magnitudes. We use the MLCS magnitudes for 
the R98 data, both including and excluding the unclassified SN 1997ck at z = 0.97 (with 50 and 49 
SNe la, respectively). The third set are the corrected/effective stretch factor magnitudes for the 54 
Fit C SNe la of P99. In all three cases we assume that the measured magnitudes are independent. 
We also assume that SNe la at high and low z are not significantly different (Drell et al. 1999, 
Riess et al. 1999, and Wang 1999 consider the possibility and consequences of evolution), and that 
intergalactic dust has a negligible effect (Aguirre 1999 considers the effects of dust). 

Our analysis of the R98 data is similar to that described in their §4, with a few minor modifi- 
cations. For the time-independent A model we compute the likelihood function L(£Iq, Q\, Hq) for 
a range of Qq spanning the interval to 2.5 in steps of 0.1, for a range of f^A spanning the interval 



statistics (see, e.g., Bloomficld Torres & Waga 1996; Jain et al. 1998; Cooray 1999; Waga & Miceli 1999; Wang et al. 
1999), with the exception of the analysis of RQ. 



— 1 to 3 in steps of 0.1, and for a range of Hq spanning the interval 50 to 80 km s 1 Mpc 1 in steps 
of 0.5 km s- 1 Mpc" 1 . 

Our analysis of the P99 data is similar to theirs, with some modifications. Specifically, we work 
with their corrected/effective magnitudes and so need only compute a three-dimensional likelihood 
function L(Qq, Q\,Mb)- Here Mb is their "Hubble-constant-free" B band absolute magnitude 
(see P99) related to Hq by Mb = —19.46 — 51og 10 i?o + 25 (determined by us from results in P99). 
For the time-independent A model we compute this likelihood function for a range of Oo spanning 
the interval to 3 in steps of 0.1, for a range of f^A spanning the interval —1.5 to 3 in steps of 0.1, 
and for a range of Mb spanning the interval —3.95 to —2.95 in steps of 0.05 (which corresponds 
to about the same interval in Hq used in our analyses of the R98 data). Note that in our P99 
time-independent A model plots we do not show the likelihood for the whole £Iq — £1\ region over 
which we have computed it. 

The spatially-flat time-variable A model we consider is the scalar field model with potential 
V{4>) oc (j)~ a , a > 0. When a — ► the model tends to the flat-constant-A case and when a — > oo 
it approaches the Einstein-de Sitter model. It is discussed in detail in PR, RP, and RQ, and we 
derive the predicted distance moduli for the SNe using expressions given in these papers. 

As shown in these papers, the scalar field behaves like a fluid with a constant (but different) 
equation of state in each epoch of the model. For instance, in the CDM and baryon dominated 
epoch of the model, it obeys the equation of state p^ = w^p^ (relating the pressure and energy 
density of the scalar field), where 

^ = -dh' (1) 

see, e.g., RQ eq. (2) or see Zlatev et al. (1999) for a more recent derivation. We shall also have 
need for an average equation of state parameter, used by Perlmutter et al. (1999b) and Wang et 
al. (1999), 

J Q a ° daVt^ajw^a) 

J dail^a) 

where fi^ is the scalar field density parameter and a is the scale factor (with ao being the present 
value) . 

In this model, for the R98 data we compute the likelihood function L(£Iq, a, Hq), and for the 
P99 data the likelihood function L(£Iq, a, Mb)- In both cases we evaluate the likelihood function 
for a range of 0,q spanning the interval 0.05 to 0.95 in steps of 0.025, for a range of a spanning the 
interval to 8 in steps of 0.5, and for the same range of Hq or Mb as in the time-independent A 
cases discussed above. 

We marginalize these three-dimensional likelihood functions by integrating over Hq (for the 
R98 data) or Mb (for the P99 data) and derive two-dimensional likelihood functions, L(Qq,Qj^) 
for the constant A model and L(ttQ,a) for the spatially-flat time-variable A case. These two- 
dimensional likelihood functions are used to derive highest posterior density limits (see Ganga et 



al. 1997 and references therein) in the (CIq, Q\) or (Qq, a) planes. In what follows we consider 1, 
2, and 3 a confidence limits which include 68.27, 95.45, and 99.73% of the area under the likelihood 
function. 

When marginalizing over a parameter or deriving a limit from the likelihood functions, we 
consider a number of different priors. We first consider a uniform prior in the parameter integrated 
over, set to zero outside the range considered for the parameter. 

We also incorporate constraints from measurements which indicate Hq = 65 ±7 km s _1 Mpc~ 
at 1 a (see, e.g., Biggs et al. 1999; Madore et al. 1999), by using the prior 



P ^ H °' v / 2^{7 km s- 1 Mpc" 1 } 6 ^ 



{H - 65 km s~ l Mpc" 1 } 2 " 
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2{7kms- 1 Mpc" 1 } 

and from measurements which indicate to = 12± 1.3 Gyr at 1 a (i.e., 0.5 Gyr added to the globular 
cluster age estimate of Chaboyer et al. 1998) 3 , by using the prior 



P{t0) = v^{L3Gyr} eXP 



{t (H ,n ,n A )-12GjrY 
2{1.3 Gyr} 2 



(4) 



with a replacing f^ in this expression in the spatially-flat time-variable A case. Note that Figure 
2 of Chaboyer et al. (1998) may be used to establish that a Gaussian prior of the form of eq. (4) is 
a good approximation to the shifted Monte Carlo globular cluster age distribution. See Ganga et 
al. (1997) for a discussion of Gaussian priors in a related context. This method of incorporating 
constraints from measurements of Hq and to is n °t identical to the classical i^o^o cosmological test 
(see, e.g., §13 of Peebles 1993). 

Finally, since f2o is a positive quantity, we also consider the non-informative prior (Berger 
1985) 4 , 

p(«o) = 7T- ( 5 ) 

Since f^o is bounded from below by the observed lower limit on the baryon density parameter 
(~ 0.01 — 0.05, depending on data used, see, e.g., Olive, Steigman, & Walker 1999), this prior does 
not result in an infinity. 



3 A more complete analysis would need to account for the uncertainty in this (0.5 Gyr) numerical value. This 
would likely weaken the effect of this prior. 

4 We thank R. Gott for emphasizing this prior; a more complete discussion of it may be found in Gott et al. (in 
preparation). I. Wasserman has noted that such a prior is probably more appropriate for a parameter that sets the 
scale for the problem, such as Ho here (see §B.2 of Drell et al. 1999; §VII of Jaynes 1968 gives a more general 
discussion) . It is however still of interest to determine how the conclusions depend on the choice of prior. 



3. Results and Discussion 

Figure 1 shows the posterior probability density distribution function (PDF) confidence con- 
tours for the time- variable A scalar field model, derived using the three data sets discussed above. 
Figure Id shows that constraints from the three different data sets are quite consistent. At 2 a a 
large region of the parameter space of these spatially-flat models (with a constant or time-variable 
A) is consistent with the SNe la data. This data favors a smaller ^o as a is increased from zero. 

A fluid with a time-independent equation of state p = wp, w < 0, has often been used to 
approximate the scalar field with potential V(<p) oc <§r a in the time-variable A model. The solid 
lines in Figure 2a show the confidence contours for such a fluid model, derived using the P99 Fit C 
data. These confidence contours are consistent with those shown in Figure 1 of Perlmutter et al. 
(1999b) (but see discussion below), Figure 10 of Wang et al. (1999), and Figure 4 of Efstathiou 
(1999). The dashed lines in Figure 2a show the exact scalar field model contours of Figure lc, 
transformed using the relation between w^ and a in the CDM and baryon dominated epoch (eq. 
[1]), which comes to an end just before the present. The two sets of contours agree near w ~ — 1, 
as they must since at w = — 1 this is the flat- A model and A does behave exactly like a fluid 
with a time-independent equation of state. However, the two sets of confidence contours differ 
significantly at larger w. Since the scalar field w^ eventually switches over to —1 (in the scalar field 
dominated epoch, RP), it is unclear what significance should be ascribed to this difference. We 
stress, however, that since the time-independent equation of state fluid model is an approximation 
to the time-variable A scalar field model, constraints on model-parameter values that are based 
on the fluid model approximation are only approximate (and possibly indicative). In passing, we 
note that Perlmutter et al. (1999b) use w e g (eq. [2]) and not w to parametrize the fluid model 
constraints. Figure 26 shows contours of constant w e s as a function of a and Slo in the time- variable 
A scalar field model. In a large part of model parameter space w e s is a sensitive function of both 
a and f^o an d hence is not the best parameter to use to describe the time-variable A scalar field 
model. 

Figure 3 shows the effects of incorporating constraints based on Hq and to measurements (eqs. 
[3, 4]). Panel o) shows that adding the Hq constraint does not significantly alter the contours 
derived from the SNe la data alone. This is expected since the value of Hq used here is very 
close to the value that is indicated by the SNe la data (see R98). However, incorporating the to 
constraint, to = 12 ± 1.3 Gyr at 1 a, does significantly shift the contours, panel b). This is because 
the SNe la data alone favor a higher t , 14.2 ± 1.7 Gyr (R98), or 14.5 ± 1.0 {0.63/ h) Gyr (P99). 

Figure 4 shows constraints on the time-variable A model, from the three different SNe la 
data sets used in conjunction with the Hq and to measurements (eqs. [3, 4]). The main effect of 
incorporating the Hq and to constraints is to increase the favored values of Qq; a weaker effect is 
the disfavoring of larger values of a. Even this extended set of data does not tightly constrain 
model-parameter values. 

Figure 5 shows the corresponding constraints on the constant A model (from the SNe la, Hq, 



and to measurements). Again, the major effect of including the Hq and to data is to increase the 
favored values of CIq. A weaker effect is that it reduces the odds against reasonable open models 
(Gorski et al. 1998), but not by a large factor. 

Figures 6 and 7 show the effects of using the non-informative prior, 1/£Iq, of eq. (5). The 
major effect is a decrease in the favored values of Qq. In the time- variable A case there is also a 
slight increase in the favored values of a (see Figure 6) while in the constant A case there is a mild 
reduction in the odds against reasonable open models (see Figure 7). If the PDF was narrower 
(i.e., if the error bars on the data were smaller), changing from the flat to the non-informative prior 
would not result in as large a change in the confidence contours. 



4. Conclusion 

Recent SNe la data do favor models with a constant or time-variable A over an open model 
without a A. However, this is not at a very high level of statistical significance. Also, the incomplete 
understanding of a number of astrophysical effects and processes (evolution, intergalactic dust, etc.) 
means that these results are preliminary and not yet definitive. 

The constraints on the time-variable A model derived here are based on the exact scalar 
field model equations of motion, not on the widely used time-independent equation of state fluid 
approximation equations of motion. 

We acknowledge the advice and assistance of R. Gott, B. Kirshner, L. Krauss, V. Periwal, 
S. Perlmutter, A. Riess, E. Sidky, M. Vogeley, and I. Wasserman, and are specially indebted to J. 
Peebles and T. Souradeep. We acknowledge helpful discussions with I. Waga who has also computed 
the PDF confidence contours for the time-variable A scalar field model, with results consistent with 
those found in this paper. We thank the referee, I. Wasserman, for a prompt and detailed report 
which helped us improve the manuscript. We acknowledge support from NSF CAREER grant 
AST-9875031. 



REFERENCES 

Aguirre, A.N. 1999, ApJ, 512, L19 

Amendola, L., Corasaniti, S., k Occhionero, F. 1999, astro-ph/9907222 

Bahcall, N.A., Ostriker, J. P., Perlmutter, S., k Steinhardt, P.J. 1999, Science, 284, 1481 

Baker, J.C., et al. 1999, MNRAS, submitted 

Banks, T., Dine, M., k Nelson, A.E. 1999, J. High Ener. Phys., 9906, 014 

Barr, S.M. 1999, Phys. Lett. B, 454, 92 

Barrow, J.D., k Magueijo, J. 1999, astro-ph/9907354 

Bartelmann, M., Huss, A., Colberg, J.M., Jenkins, A., k Pearce, F.R. 1998, A&A, 330, 1 

Bento, M.C., k Bertolami, O. 1999, gr-qc/9905075 

Berger, J.O. 1985, Statistical Decision Theory and Bayesian Analysis (New York: Springer- Verlag) , 

82 

Biggs, A.D., Browne, I.W.A., Helbig, P., Koopmans, L.V.E., Wilkinson, P.N., k Perley, R.A. 1999, 
MNRAS, 304, 349 

Binetruy, P. 1999, Phys. Rev. D, 60, 063502 

Bloomfield Torres, L.F., k Waga, I. 1996, MNRAS, 279, 712 

Bond, J.R., Jaffe, A.H., k Knox, L. 1998, Phys. Rev. D, 57, 2117 

Brax, P., k Martin, J. 1999, Phys. Lett. B, in press 

Bucher, M., k Spergel, D. 1999, Phys. Rev. D, 60, 043505 

Caldwell, R.R., Dave, R., k Steinhardt, P.J. 1998, Phys. Rev. Lett., 80, 1582 

Carroll, S.M. 1998, Phys. Rev. Lett., 81, 3067 

Chaboyer, B., Demarque, P., Kernan, P.J., k Krauss, L.M. 1998, ApJ, 494, 96 

Chiba, T., Sugiyama, N., k Nakamura, T. 1997, MNRAS, 289, L5 

Choi, K. 1999, hep-ph/9902292 

Cole, S., Weinberg, D.H., Frenk, C.S., k Ratra, B. 1997, MNRAS, 289, 37 

Colley, W.N., Gott, J.R., Weinberg, D.H., Park, C, k Berlind, A.A. 2000, ApJ, 529, in press 

Copeland, E.J., Liddle, A.R., k Wands, D. 1998, Phys. Rev. D, 57, 4686 

Cooray, A.R. 1999, A&A, 341, 653 

Donahue, M. k Voit, CM. 1999, ApJ, 523, L137 

Doroshkevich, A.G., Miiller, V., Retzlaff, J., k Turchaninov, V. 1999, MNRAS, 306, 575 

Drell, P.S., Loredo, T.J., k Wasserman, I. 1999, ApJ, in press 

Efstathiou, G. 1999, MNRAS, submitted 



10 



Efstathiou, G., Sutherland, W.J., & Maddox, S.J. 1990, Nature, 348, 705 

Falco, E.E., Kochanek, C.S., & Muhoz, J.A. 1998, ApJ, 494, 47 

Ferreira, P.G., & Joyce, M. 1998, Phys. Rev. D, 58, 023503 

Fischler, W., Ratra, B., & Susskind, L. 1985, Nucl. Phys. B, 259, 730 

Freese, K., Adams, F.C., Frieman, J.A., k Mottola, E. 1987, Nucl. Phys. B, 287, 797 

Freudling, W., et al. 1999, ApJ, 523, 1 

Frieman, J.A., Hill, C.T., Stebbins, A., & Waga, I. 1995, Phys. Rev. Lett., 75, 2077 

Frieman, J.A., & Waga, I. 1998, Phys. Rev. D, 57, 4642 

Ganga, K., Ratra, B., Gundersen, J.O., & Sugiyama, N. 1997, ApJ, 484, 7 

Ganga, K., Ratra, B., Lim, M.A., Sugiyama, N., & Tanaka, S.T. 1998, ApJS, 114, 165 

Ganga, K., Ratra, B., & Sugiyama, N. 1996, ApJ, 461, L61 

Garcia-Berro, E., Gaztahaga, E., Isern, J., Benvenuto, O., & Althaus, L. 1999, astroph/9907440 

Garnavich, P.M., et al. 1998, ApJ, 509, 74 

Garriga, J., Livio, M., & Vilenkin, A. 1999, astro-ph/9906210 

Giovannini, M. 1999, astro-ph/9903004 

Gorski, K.M., Ratra, B., Stompor, R., Sugiyama, N., &: Banday, A.J. 1998, ApJS, 114, 1 

Gorski, K.M., Ratra, B., Sugiyama, N., & Banday, A.J. 1995, ApJ, 444, L65 

Gott, J.R. 1982, Nature, 295, 304 

Gott, J.R. 1997, in Critical Dialogues in Cosmology, ed. N. Turok (Singapore: World Scientific), 
519 

Gratton, R.G., Fusi Pecci, F., Carretta, E., Clementini, C, Corsi, C.E., & Lattanzi, M. 1997, ApJ, 
491, 749 

Hu, W., Eisenstein, D.J., Tegmark, M, & White, M. 1999, Phys. Rev. D, 59, 023512 

Hoflich, P., Wheeler, J.C., & Thielemann, F.K. 1998, ApJ, 495, 617 

Jain, D., Panchapakesan, N., Mahajan, S., & Bhatia, V.B. 1998, ApJ, submitted 

Jaynes, E.T. 1968, IEEE Trans. Sys. Sci. Cyb., SSC-4, 227 

Kamionkowski, M. & Kosowsky, A. 1999, Ann. Rev. Nucl. Part. Sci, in press 

Kim, J.E. 1999, J. High. Ener. Phys, 9905, 022 

Kolda, C. & Lyth, D.H. 1999, Phys. Lett. B, 458, 197 

Kravtsov, A.V., & Klypin, A.A. 1999, ApJ, 520, 437 

Liddle, A.R., & Scherrer, R.J. 1999, Phys. Rev. D, 59, 023509 

Lineweaver, C.H., & Barbosa, D. 1998, ApJ, 496, 624 



11 



Lucchin, F., k Matarrese, S. 1985, Phys. Rev. D, 32, 1316 

Madore, B.F., et al. 1999, ApJ, 515, 29 

Masiero, A., Pietroni, M., k Rosati F. 1999, Phys. Rev. D, submitted 

Meneghetti, M., Bolzonella, M., Bartelmann, M., Moscardini, L., k Tormen, G. 1999, MNRAS, 
submitted 

Olive, K.A., Steigman, G., k Walker, T.P. 1999, Phys. Rep., in press 

Overduin, J.M. 1999, ApJ, 517, LI 

Ozer, M. 1999, ApJ, 520, 45 

Ozer, M., k Taha, M.O. 1986, Phys. Lett. B, 171, 363 

Peebles, P.J.E. 1984, ApJ, 284, 439 

Peebles, P.J.E. 1993, Principles of Physical Cosmology (Princeton: Princeton University Press) 

Peebles, P.J.E., k Ratra, B. 1988, ApJ, 325, L17 (PR) 

Peebles, P.J.E., k Vilenkin, A. 1999, Phys. Rev. D, 59, 063505 

Periwal, V. 1999, astro-ph/9906253 

Perlmutter, S., et al. 1997, ApJ, 483, 565 

Perlmutter, S., et al. 1999a, ApJ, 517, 565 (P99) 

Perlmutter, S., Turner, M.S., k White, M. 1999b, Phys. Rev. Lett., 83, 670 

Perrotta, F., k Baccigalupi, C. 1999, Phys. Rev. D, 59, 123508 

Ratra, B. 1991, Phys. Rev. D, 43, 3802 

Ratra, B. 1992, Phys. Rev. D, 45, 1913 

Ratra, B., Ganga, K., Stompor, R., Sugiyama, N., de Bernardis, P., k Gorski, K.M. 1999a, ApJ, 
510, 11 

Ratra, B., k Peebles, P.J.E. 1988, Phys. Rev. D, 37, 3406 (RP) 

Ratra, B., k Peebles, P.J.E. 1994, ApJ, 432, L5 

Ratra, B., k Peebles, P.J.E. 1995, Phys. Rev. D, 52, 1837 

Ratra, B., k Quillen, A. 1992, MNRAS, 259, 738 (RQ) 

Ratra, B., Stompor, R., Ganga, K., Rocha, G., Sugiyama, N., k Gorski, K.M. 1999b, ApJ, 517, 
549 

Reid, I.N. 1997, AJ, 114, 161 

Riess, A.G., et al. 1998, AJ, 116, 1009 (R98) 

Riess, A.G., Filippenko, A.V., Li, W., k Schmidt, B.P. 1999, AJ, in press 

Rocha, G. 1999, in Dark Matter in Astrophysics and Particle Physics 1998, ed. H.V. Klapdor- 
Kleingrothaus k L. Baudis (Bristol: Institute of Physics Publishing), 238 



12 



Rocha, G., Stompor, R., Ganga, K., Ratra, B., Piatt, S.R., Sugiyama, N., & Gorski, K.M. 1999, 
ApJ, 525, in press, astro-ph/9905127 

Sahni, V., & Starobinsky, A. 1999, Int. J. Mod. Phys. D, in press 

Santiago, D.I., & Silbergleit, A.S. 1998, Phys. Lett. A, submitted 

Stompor, R., Gorski, K.M., & Banday, A.J. 1995, MNRAS, 277, 1225 

Suntzeff, N.B., et al. 1999, AJ, 117, 1175 

Turner, M.S., & White, M. 1997, Phys. Rev. D, 56, 4439 

Umeda, H., Nomoto, K., Kobayashi, C, Hachisu, I., k Kato, M. 1999, ApJ, 522, L43 

Uzan, J.-P. 1999, Phys. Rev. D, 59, 123510 

Waga, I., & Miceli, A.P.M.R. 1999, Phys. Rev. D, 59, 103507 

Wang, L., Caldwell, R.R., Ostriker, J.P., & Steinhardt, P.J. 1999, ApJ, submitted 

Wang, Y. 1999, ApJ, submitted 

Wetterich, C. 1995, A&A, 301, 321 

Yamamoto, K., & Bunn, E.F. 1996, ApJ, 464, 8 

Zlatev, I., Wang, L., & Steinhardt, P.J. 1999, Phys. Rev. Lett., 82, 896 



This preprint was prepared with the AAS IATpjX macros v5.0. 



13 



Figure Captions 

Fig. 1— PDF confidence contours for the spatially-flat time- variable A scalar field model, with 
potential V((p) oc 4>~ a , derived using the three SNe la data sets. The a = axis corresponds to the 
spatially-flat time-independent A case. Confidence contours in panels a) — c) run from —2 to +3 a 
starting from the lower left hand corner of each panel. Panel a) shows those derived from all the 
R98 SNe, while b) is for R98 SNe excluding the z = 0.97 one, and c) is for the P99 Fit C data set. 
Panel d) compares the ±2 a limits from the three data sets: all R98 SNe (long-dashed lines); R98 
SNe excluding the z = 0.97 one (short-dashed lines); and the P99 Fit C SNe (dotted lines). 

Fig. 2.- a) PDF confidence contours derived from the P99 Fit C SNe, for a spatially-flat fluid 
model with equation of state p = wp (solid lines show contours from —3 to +3 a starting from the 
lower left hand corner) and for the spatially-flat time-variable A scalar field model with potential 
V(4>) oc (j)~ a which behaves like a fluid model with equation of state ps = w^p^ in the CDM and 
baryon dominated epoch (dotted lines show contours from —2 to +3 a starting from the lower 
left hand corner). The scalar field model contours were derived by using eq. (1) to transform 
the vertical axis of Figure lc. The fluid model likelihood function was computed for the full 
range of w and Qq shown. The dot-dashed lines bound the w^ — 0,q region that corresponds to 
the a — Qq region of Figure lc over which the scalar field model likelihood was computed, b) 
Contours of constant tu c g (eq. [2]) in the a — 0,q plane for the time- variable A scalar field model 
with potential V((f>) oc <p~ a . Starting at the lower left hand corner, the contours correspond to 
w eS = -0.9, -0.8, -0.7, -0.6, -0.5, -0.4, -0.35, -0.3, and -0.25. 

Fig. 3.- PDF confidence contours derived from the P99 Fit C SNe, for the spatially-flat time- 
variable A scalar field model with potential V(4>) oc <ft~ a . The dotted lines in panels a) — c) show 
the —2 to +3 a contours of Figure lc. The solid lines in these panels show the corresponding 
contours derived from the SNe data in conjunction with: a) Hq measurements (using eq. [3]); b) 
to measurements (using eq. [4]); and c) Ho and to measurements (using eqs. [3, 4]). Panel d) 
compares the ±2 a confidence contours from the SNe data in conjunction with: the Ho constraint 
(long-dashed lines); the to constraint (short-dashed lines); and both constraints (dotted lines). 

Fig. 4.- PDF confidence contours for the spatially-flat time-variable A scalar field model, with 
potential V((p) oc (f>~ a , derived using the three SNe la data sets. Panel a) shows those derived from 
all the R98 SNe, while b) is from the R98 SNe excluding the z = 0.97 one, and c) is for the P99 
Fit C data. Solid (dotted) lines in a) — c) are the 1, 2, and 3 a contours from the SNe data with 
(without) the Ho and to constraints (eqs. [3, 4]); the dotted lines here are the solid lines in Figure 
la — c. Panel d) compares the ±2 a limits from the Ho and to constraints used in conjunction with: 
all R98 SNe (long-dashed lines); R98 SNe excluding the z = 0.97 one (short-dashed lines); and the 
P99 Fit C SNe (dotted lines). 

Fig. 5.- PDF confidence contours for the time-independent A model, derived using the three SNe 
la data sets. Panel a) shows those derived from all the R98 SNe, while b) is from the R98 SNe 
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excluding the z = 0.97 one, and c) is for the P99 Fit C SNe. Solid (dotted) lines in a) — c) are 
the 1, 2, and 3 a contours from the SNe data with (without) the Hq and to constraints (eqs. [3, 
4]). Panel d) compares the ±2 a limits from the Hq and to constraints used in conjunction with: 
all R98 SNe (long-dashed lines); R98 SNe excluding the z = 0.97 one (short-dashed lines); and the 
P99 Fit C SNe (dotted lines). In all panels, models with parameter values in the upper left hand 
corner region bounded by the diagonal dot-dashed curve do not have a big bang. The horizontal 
dot-dashed line demarcates models with a zero A and the diagonal dot-dashed line running from 
f&A = 1 to r^o = 2 corresponds to spatially-flat models. 

Fig. 6.- PDF confidence contours (1, 2, and 3 a) derived from the P99 Fit C SNe, for the spatially- 
flat time-variable A scalar field model. Panel a) (b)) ignores (accounts for) the Hq and to constraints 
(eqs. [3, 4]). The solid (dotted) lines use the non-informative 1/fio (flat) prior. The dotted lines 
in a)(b)) are the same as the solid lines in Figure lc (4c). 

Fig. 7.- PDF confidence contours (1, 2, and 3 a) derived from the R98 SNe excluding the z = 0.97 
one, for the time- independent A model. Panel a) (&)) ignores (accounts for) the Hq and to constraints 
(eqs. [3, 4]). The solid (dotted) lines use the non-informative 1/CIq (flat) prior. The dotted lines 
in a) (6)) are the same as the dotted (solid) lines in Figure 56. The dot-dashed lines are described 
in the caption of Figure 5. In the non-informative prior cases the likelihood function is computed 
down to f^o = 0.01. 
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